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SUMMARY 

A brief study was made to assess the applicability of the Newton- Raphson ^gital 
comDuter program as a routine technique for extracting aerodynamic derivatives from 
S tests of^ liftTng body types of vehicles. Lateral-directional flight data from flight 
tests of th ' HL-tO lifting body research vehicle were utilized. The results, in gener , 
show the computer program to be a reliable and expedient means for extracting deriva- 
tives for this class of vehicle as a standard procedure. Tms result was true even wh 
Sflity augmentation was used. An - a priori- weighting option available in the com- 
puter program was found to be a desirable feature. As a resuU of the study. a credib 
set of HL-10 lateral-directional derivatives was obtained from flight data. These deriv 
atives are compared with results from wind-tunnel tests. 

INTRODUCTION 


Prior to the start of the M2-F3 lifting body research vehicle flight program (ref. I ), 
the need existed for assessing the reliability of using the Newton -Raphson digital com- 
a. , i#v rpfc o a \ oc the nrinoipal means o r derivative extraction tor future 

ft A wa. therefore conducted tn wh.eh prevt^l, 

obtained lateral-directional flight data of the HL-10 lifting body were used. This vehicle 
Zs chosen because 1. was well documented in terms of wlnd-tunne “ > re '- « ^ 

earlier flight evaluations of its stability and control charac teristics (ref. 6). In additi , 
time histories of transient maneuvers were readily available in digitized for . 

Two ground rules were used to conduct this study. First, various options available 
in the computer program were held fixed, and, second, once these options^ were estab- 
lished each maneuver was run one time only. These rules were established to ^ 
tain whether an acce. able set of flight derivatives could be obtained by using only one 
computer run per maneuver. 

The resulting set of HL-10 lateral-directional derivatives obtained is presented nnd 


compared with wind-tunnel data of reference 5. In addition, the effects of using the 
computer program a priori option are included. 
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SYMBOLS 


Derivatives are presented as standard NASA coefficients of forces and moments. A 
right-hand sign convention is used to determine the direction of all forces, moments, 
angular displacements, and velocity. 

Physical quantities are given in both the International System of Units (SI) and paren- 
thetically in U.S. Customary Units. All measurements were taken m U.S. Customary 
Units. Conversion factors can be found in reference 7. 

A stability matrix, P x P 

a transverse acceleration, g 

1 y 

3 control matrix, P * Q 

b reference body span, m (ft) 

C transformation matrix, P x P 

F force, N (lb) 

q partition of matrix relating the state vector to the observation 

vector, (R - P) x P 

2 2 

g acceleration due to gravity, 9. x m/sec (32.2 ft/ sec ) 

3 partition of matrix relating the control vector to the observation 

vector, <R - P) x Q 


l XZ 


M 

T\ 


m 

O 


identity matrix 

2 2 

rolling moment oi inertia, kg-m (slug-ft > 

2 2 

product of inertia, kg-m (slug-ft ) 

2 2 

yawing moment of inertia, kg-m (slug-ft ) 

Mach number 
moment, m-N (ft-lht 
mass, kg (slugs) 
null matrix 
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■■'V 


li 


p 

p 

Q 

q 

R 

r 

S 

u 

Y 
X 
x 

Y 

Y 
Z 
a 

P 
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a 


number of state variables 

rolling rate, rad/sec or deg/sec 

number of control variables 

2 2 

dynamic pressure, N/m (lb/ft ) 

number of observation variables 

yawing rate, rad/sec or deg/see 

2 . 2 

reference planform area, m (it ) 
control vector (Q x 1) 
velocity, m/sec (ft/ sec) 

X-axis 

state vector, P x 1 
Y-axis 

observation vector, R x 1 
Z-axis 

angle of attack, deg 
angle of sideslip, deg 
aileron deflection, deg 

rudder deflection, deg 

constant control deflection, rad or deg 


<P 

C, 


pitching attitude, deg 
angle of bank, deg 


M 


X 


rolling-moment coefficient, qg^ 




yawing -moment coefficient, ^g^ 

t y 

side-force coefficient, trg- 
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L §^c 
<* ‘x h 


N . - 

P l z j3 


Y = — C 
P mV i 


(3 


qSb 2 c 

L p 2 VI C l 
y X P 


N - S§l 2 _- c 
P 2 VI z C n p 



L fi fi = ^ C l 

6 a ,6 r J X l 6 ,6 
a’ r 


isi = 

^ 6^,6 \ry n* 


a' r A Z 6 a ,6 r 


Y = - 3 §_ C 

6 a’ 6 r mV Y 6 a ,6 r 


L =^-C 
r 2 VI x C / r 


V - q Sb r 

N r-5VT C n r 


Y = sin a 
P 

Subscripts: 


X 

Y 

Z 

6 a‘ 6 r , 6 o ’ P ’ r,/3 


X-axis component 
Y-axis component 
Z-axis component 

partial differentiation with respect to subscripted variable 


A dot over a symbol signifies a derivative with respect to time. 


METHOD OF ANALYSIS 


A digital computer program was used to identify lateral-directional sets of deriva- 
f . A ftS fUahTdata This program uses a modification of the Newcon-Raphson method 
'" d calwTteNa^fon-Raph..* program. Detailed ^cussion of the .pro- 
theory and applications, is given in references 2 to 4. The appendix 
contains'the^set of equations fmode.) used to identify the derivatives for this report. 

The method is « —e J "Ss 

alternatives to minimize the difference between computed and measured histoi ics. 

The output time histories are assumed to contain noise, but the input (coptro ) 
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histories are defined as being noise free. The input time histories used were those tor 
the ailerons (differential elevens) and rudder. The output time histories were those lor 
roil rate, yaw rate, sideslip angle, bank angle, and lateral acceleration. 

An option used, called "a priori" (ref. 2), allowed the startmg set of derivatives to 
be weighed, which tended to hold the resulting derivatives near their start mg value if no 
information about them was contained in the maneuver. The procedure used is to mc ude 
a penalty for departure from the assumed starting value. For this study, wind-tunne 
predictions were used as starting values. 

DESCRIPTION OF VEHICLE 

The HL-10 is a research vehicle of the lifting body class. Its basic aerodynamic 
shai is that a thick, 74 sweep, delta plahtom with ne^ive camber ahd Uirec ■ aH _ 
vertical fins. Pertinent physical characteristics are presented in table I, and a three 
view drawing is shown in figure 1. 

The primary controls consisted of elevens and a rudder. The elevens provide-d both 
pitch and roll control, damping augmentation, and trim control. The judder located 
on the center vertical fin, was split so as to lie operated botn as a ruddei and a sp 
brake. Secondary control surfaces were located on the inboard and outboard trailing 
edees of the tip fins and upper surfaces of the elevons. These su daces serve< 
position flaps and were deployed either open for transonic flight or e ose< or su sonic 

flight. 

The primary control surfaces were actuated by irreversible hydraulic systems and 
accepted 3 commands from both the pilot and the stability augmentation system (SAS). The 
SAS consisted of a three-axis rate feedback system with pilot adjustable gams. 

Further vehicle description can be found in references (>, *, 9, and 10. 

INST RU MEN T ATION 

Conventional NASA instrumentation was used to determine all flight quantities of 
interest (ref. fi). Data were acquired by means of a 9-bit pulse code mot u a ion y 
and were telemetered in real time to a ground station. Measurement aceuiacies «»f 
these quantities are detailed in reference <>. 


FLIGHT TESTS 
P rocedures 

Frequent weight and balance measurements were made to verify the 
vehicle center of gravity. Moments of inertia were determined experimentally before 
the initial flight by means of an inertia swing (ref. 11). The inertia estimate was 
updated analytically whenever the mass distribution changed. 

Like other lifting bodies, the HL-10 vehicle was air-launched from a modified B-52 


airplane at an altitude of approximately 14,000 meters (45,000 feet) and a Mach number 

of O 67? After launch, the pilot flew a preplanned flight profile. The unpowered, or 

glide flights lasted less than 4 minutes and were usually made below a Mach number o 
0 7 ’For powered flights, the engine was lit immediately after launch, angle oi attack 
waVinc^eied S pdn altfeide. and then the vehicle was pushed over to increase Mach 
^mS ?he peered portion of the flight, which usually lasted from 90 seconds to 180 
•seconds was made in the transonic configuration. A change to the subsonic configura 
tion was’ made when the Mach number decreased to about 0.7. The altitude at this , 

was about 9150 meters (30, 000 feet). Most of the stability and control data were obk 

after engine burnout. 

In general, maneuvers from which data were obtained were performed at altitudes 
above approximately 6100 meters (20,000 feet) to provide the pilot with enough time to 
set up for the final approach and landing. The trajectories flown precluded steady 
fli^at conditions. In addition, the SAS was generally used throughout the flight profile. 

MANEUVERS 

The maneuvers used had previously been performed for the purpose of analog 
derivative extraction (ref. 6). These were the standard doublet types of maneuvers 
which have been used (ref. 8) to excite the airplane’s transient characteristics. The 
SAS was generally used to insure satisfactory handling qualities during maneuvers. 
However, damper gains were usually at or less than 0.5 deg/deg/sec. All maneuvers 
used in this study were performed with the vehicle in the transonic configuration. Con 
trol derivatives obtained from a maneuver were used only if the maneuver containe a 
pilot input from that control. 

RESULTS AND DISCUSSION 


Summary of Results 

For this study an attempt was made to analyze 42 lateral-directional maneuvers 
(i.e., all that were available in convenient digitized form). Acceptable sets of deriva- 
tives were obtained in 26 cases. Of the maneuvers which yielded unacceptable results, 

13 converged to a set of derivatives with an unacceptable fit error, while three were 
not convergent. Of the 13 convergent cases, two cont. \ed atmospheric turbulence, 
seven were very weak in amplitude, and the remaining rour were of very large amplitude. 
It should be noted that four of the unacceptable cases were later found to give acceptable 
results after changes were made in the program options. To follow the ground rules 
and keep the study consistent, these four cases are not presented. 

Figures 2, 3, and 4 present time histories of the typical, best, and worst (but still 
acceptable) matches obtained. 
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Comparison of Flight and Wind-Tunnel Data 
Figures 5 to 9 compare the flight derivatives with wind-tunnel predictions. All 


derivatives are plotted as functions of angle of attack at selected wind-tunnel test Mach 
numbers. Wind-tunnel values are based on data from reference •>. 

An overview of figures 5 to 9 indicates that the effective dihedral derivative, , 

and the directional stability derivative, C , agree well with wind-tunnel data. The 

li 

side-force derivative, C y , is often lower than wind-tunnel predictions. Aileron 

control derivatives are approximately at the same level \s wind-tunnel data. However, 

0 f i 9 ifior is more proverse than wind-tunne. piedietions at 

n 6 n 

high angles of attack. The side force due to aileron, C y ^ , from flight was not 

cl 

determined (i.e., held fixed) in this study. The rudder effectiveness, , substan- 

r 

tinted wind-tunnel predictions. The rolling moment due to rudder, , was approxi- 

r 

matelv at the same level as wind-tunnel data, although it tended to have a good deal of 
"c atteV: indicattag that it was relatively hard to identify. The side force due to rudder, 
C , was consistently lower than the wind-tunnel data* 

Y o 

r 

Damping derivatives are generally more difficult to identify than other derivatives 
because of their small magnitude and the use of the SAS during transient maneuvers 
The damping derivatives (C , C , C , and C n ) show enough scatter to suppoit 

c p c r p r 

this. 


Effect of the A Priori Option 

Figure 10 compares the derivatives obtained from High! data with and without the 
a priori option for a Mach number of 1.2. The general level ot the major panamete^ 

J r 1 C C C and C ) is not significantly affected In asing ti.c option. 




a 


Scatter in the damping derivatives (C ; , C 


n 

P P 


and C ) has been reduced, 


although there is too much scatter in these derivatives obtained without a priori to deter- 
mine S those with a priori were biased. The weighting on C y appears to he too high, 

suc'h that this derivative was biased toward the wind-tunnel value. 

It has been found that the option tends to help C onverg«*nce--especia lly when the 
maneuver does not contain sufficient information to define all the derivatives well. 1 
option a as most successfully used to obtain derivatives from maneuvers with »tahilit\ 
augmentation anil maneuvers with weak pilot control inputs. 


( 


CONCLUDING REMARKS 





A study was made of the feasibility of using the Newton-Raphson digital computer 
program to identify lifting body derivatives from flight data. The study showed that e 
Newton-Raphson program was a reliable and expedient means of extracting these deriv- 
atives on a routine basis. This was true even when stability augmentation was used. 
The a priori option w ? as shown to be a desirable feature. 

In general, a satisfactory set of HL-10 lifting body lateral-directional derivatives 
was obtained from the flight data. 


Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif., November 2, 1973 
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APPENDIX 


EQUATIONS OF MOTION MECHANIZED IN THE NEWTON-RAPHSON 
DIGITAL COMPUTER PROGRAM 

The following state equations were used in the basic model for this study. 

Cx = Ax + Bu 

X=[-i]5 + C-g-]H 


where x, x, u , and £ are Mme varying. 
For the lateral-directional mechanization. 
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TABLE 1.- PHYSICAL CHARACTERISTICS OF THE HL-10 LIFTING BODY VEHICLE 


(a) Reference areas and lengths 



Body - 22 

Reierence planform area, m (ft ) 

Length, m (ft) 

Span, m (ft) 

2 

Aspect ratio, b /S 

Elevons (two) - 

Area, each, m 2 (ft ) 

Span, each, m (ft) 

Chord: 

Root, m (ft) 

Tip, m (ft) • • 

Elevon flap (two) - 

2 2 

Area, each, m (ft ) 

Span, each, m (ft) 

Chord: 

Root, m (ft) • • 

Tip, m (ft) 

Vertical stabilizer - 

Area, m 2 (ft 2 ) 

Height, m (ft) 

Chord: 

Root, m (ft) 

Tip, m (ft) 

Leading-edge sweep, deg 

Rudders (two) - 

2 2 

Area, each, m (ft ) 

Height, each, m (ft) 

Chord, m (ft) 

Outboard tip fin flaps (two) - 
2 2 

Area, each, m (ft ) 

Height at hinge line, m (ft) 

Chord perpendicular to hinge line, m (ft) , 
Inboard tip fin flaps (two) - 

Area, m 2 (ft 2 ) 

Height at hinge line, m (ft) 

Chord perpendicular to hinge line, m (ft) . 


14 . 9 ( 160 ) 
6.45 ( 21 . 17 ) 
4 . 15 ( 13 . 60 ) 

1 . 156 


1 . 00 ( 10 . 72 ) 
1. 09 ( 3 . 58 ) 

0 . 59 ( 1 . 93 ) 
1.24 ( 4 . 06 ) 


0.70 ( 7 . 50 ) 
1.09 ( 3 . 58 ) 

0.48 ( 1 . 58 ) 
0 . 80 ( 2 . 63 ) 


1 . 47 ( 15 . 80 ) 
1. 53 ( 5 . 02 ) 

1.32 ( 4 . 32 ) 
0.60 ( 1 . 97 ) 
25 


0.41 ( 4 . 45 ) 
1.26 ; 4 . 12 ) 
0 . 33 ( 1 . 08 ) 


0 . 35 ( 3 . 77 ) 
1. 37 ( 4 . 50 ! 
0 . 76 ( 2. 48 


0 . 23 ( 2 . 48 ) 
1.01 ( 3 . 31 ) 
0.23 ( 0 . 75 ) 
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Figure *. Comparison of derivatives obtained from flight data with wind-tunnel 
predictions for a Mach number of 1.2. 
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Figure 10. Comparison of derivatives obtained from flight data with and without 
the a priori option for a Mach number of 1.2. 
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